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Abstract 

Motivated by the discovery of errors in six of the 135 diagrams in the published 
five-loop expansions of the /3-function and the anomalous dimensions of the 
0(n)-symmetric </> 4 -theory in D = 4 — e dimensions we present the results of 
a full analytic reevaluation of all diagrams. The divergences are removed by 
minimal subtraction and e-expansions are given for the critical exponents rj, u, 
and oj up to order e 5 . 
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1) During the last two decades, much effort has been invested into studying the 
scalar quantum field theory with 4 -interaction. On the one hand, such a theory 
describes correctly many experimentally observable features of critical phenomena. 
Field theoretic renormalization group techniques in D = 4 — e dimensions 0, [| |J 
combined with Borel resummation methods of the resulting e-expansions || led to 
extremely accurate determinations of the critical exponents of all 0(n) universality 
classes. The latter requires the knowledge of the asymptotic behaviour of perturbation 
series in four dimensions which is completely known in this theory |J . Apart from such 
important applications, the 4 -theory, being the simplest renormalizable quantum 
field theory in the four dimensional space-time, has been an ideal ground for testing 
new methods of calculating Feynman diagrams and for studying the structure of 
perturbation theory. 

The RG functions of the (/> 4 -theory were first calculated analytically in four di- 
mensions using dimensional regularization J7| and the minimal subtraction (MS) 
scheme || in the three- and four-loop approximations in Refs. || and [fT0| . The 
critical exponents were obtained as e-expansions || up to terms of order e 3 and e 4 . 

The five-loop anomalous dimension of the field and the associated critical expo- 



nent 7] to order e were determined analytically in [II]. The five- loop /3-function and 



the anomalous dimension of the mass were given in Ref. [12]. However, three of the 
124 four-point diagrams contributing to the /^-function at the five-loop level could 
be evaluated only numerically. The analytic calculation of the /5-function was finally 
completed in ||13|| . The ensuing e-expansions for the critical exponents were obtained 
up to order e 5 in ||14j| . 

Intending further applications, the Berlin group of the authors undertook an in- 



dependent recalculation of the perturbation series of Refs. [[□], |I2]| , using the same 
techniques, and discovered errors in six of the 135 diagrams. This meant that the 



subsequent results of [^, [14| were also incorrect. When visiting the Moscow group 
the errors were confirmed and we can now jointly report all expansions in the correct 
form. 

2) We consider the 0(n)-symmetric theory of n real scalar fields (f) a (a = 1, 2, . . . , n) 
with the Lagrangian 

in an euclidean space with D = 4 — e dimensions. The bare (unrenormalized) coupling 
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constant g-Q and mass m-Q are expressed via renormalized ones as 



g B = H e Z g g = // 



Z m 27Jl 



-m 



(2) 



Here [i is the unit of mass in dimensional regularization and Z4, Z2, Z m 2 are the 
renormalization constants of the vertex function, propagator and mass, respectively, 
with Z$2 being the renormalization constant of the two-point function obtained from 
the propagator by the insertion, in all possible ways, of the vertex (— 2 ) |TIJ. In the 
MS-scheme the renormalization constants do not depend on dimensional parameters 
and are expressible as series in 1/e with purely (^-dependent coefficients: 
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E 

fc=i 



ZiM 



(3) 



The /3-function and the anomalous dimensions entering the RG equations are ex- 
pressed in the standard way as follows: 



e d g 

2 9+ Ih^ 2 



d mm 



In 



We also use the relations 



d ln/x 
d In Z. 



d In fx 2 



SB 
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d In Z m 2 
d In fi 2 
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2 d g 
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2,4,. 



(4) 

(5) 
(6) 



Pis) = 9[ 2 i2{g) - 14(g)} , jm(g) = 12(g) - i<p(g), 



(7) 



which follow from the relations between renormalization constants implied by (0) and 
are useful for the calculations of (3(g) and i m (g). 

To determine all RG functions up to five loops we calculate the five-loop ap- 
proximation to the three constants Z 2 , Z4 and Z^. The constant Z 2 contains the 
counterterms of the 11 five- loop propagator diagrams. Two of them were calculated 



erroneously in Ref. [11]. The constant Z^ receives contributions from 124 vertex dia- 
grams. Of these diagrams, 90 contribute to Z^p, after appropriate changes of combina- 
torial factors. Four of the 124 counterterms were calculated erroneously in Ref. |T2 



In the present paper we have used the same methods as in the previous 
works |ll|, to calculate the counterterms from the dimensionally regularized Feyn- 
man integrals, namely, the method of infrared rearrangement |15[ , the Gegenbauer 
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polynomial x-space technique (GPXT) [|T(J, the integration- by-parts algorithm ||17|| , 
and the .Reoperation Three diagrams were calculated analytically first in [13j 
by using the so-called method of uniqueness, later the same results were obtained 
for them by using the Gegenbauer polynomials in x-space together with several non- 
trivial tricks [|19|]. A detailed description of the calculations including the diagramwise 
results will be presented in a separate publication. 

The analytic results of our recalculation of the five-loop approximations to the 
RG functions /3(g), 72(g) and 7 m (g) are [C( n ) is the Riemann ^-function]: 

13(g) = £[n + 8] - £[3n+14] 

[33n 2 + 922n + 2960 + C(3) ■ 96(5n + 22)] 
[-5n 3 + 6320n 2 + 80456n + 196648 
+ C(3) • 96(63n 2 + 764n + 2332) 
-C(4) -288(5n + 22)(n + 8) 
+ C(5) • 1920(2n 2 + 55n + 186)] 
+T24il6 t 13n4 + 12578n3 + 808496n 2 + 6646336n + 13177344 

+ C(3) • 16(-9n 4 + 1248n 3 + 67640n 2 + 552280n + 1314336) 
+ C 2 (3) • 768(-6n 3 - 59n 2 + 446n + 3264) 

- C(4) • 288(63n 3 + 1388n 2 + 9532n + 21120) 

+ C(5) • 256(305n 3 + 7466n 2 + 66986n + 165084) 

- C(6)(n + 8) • 9600(2n 2 + 55n + 186) 

+ C(7) • 112896(14n 2 + 189n + 526)] , (8) 



72(5) = 9 i(n + 2) - |l(n + 2)[n + 8] 



+ 5TM( n + 2 ) [5(-n 2 + 18n + 100)] 



36 v*-r*/ 432 
4 

5184 

5 

186624 

3 P-rn 2 



3—'n + 2)[39n 3 + 296n 2 + 22752n + 77056 



- C(3) • 48(n 3 - 6n 2 + 64n + 184) 
+ C(4) • 1152(5n + 22)] , (9) 



7m (<7) = |(n + 2) - j£(n + 2)[5] + fj(n + 2)[5n + 37] 



36 V" ' ^/L"J ' 72 

9 4 
15552 

,2 



3—(n + 2) [-n 2 + 7578n + 31060 



+ C(3) -48(3n 2 + 10n + 68) 
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+ C(4) • 288(5n + 22)] 
+ 37§J 3 g(n + 2)[21n 3 + 45254n 2 + 1077120n + 3166528 

+ C(3) • 48(17n 3 + 940n 2 + 8208n + 31848) 
- C 2 (3) • 768(2n 2 + 145n + 582) 
+ C(4) • 288(-3n 3 + 29n 2 + 816n + 2668) 
+ C(5) • 768(-5n 2 + 14n + 72) 

+ C(6) • 9600(2n 2 + 55n + 186)] . (10) 
For n = 1 the series have the numerical form: 

P( g ) = 1.5 g 2 - 2.833 g 3 + 16.27 g A - 135.8 g 5 + 1424.2841 g 6 , (11) 

72 = 0.0833 g 2 - 0.0625 g 3 + 0.3385 g A - 1.9256 g 5 , (12) 

7m = 0.5 g - 0.4167 g 2 + 1.75 g 3 - 9.978 g 4 + 75.3778 g 5 . (13) 

Note that the five-loop coefficients have changed by about 0.3 % for the /3-function, 
by about 9 % for j m , and by a factor of three for 72 in comparison with the wrong 



results of Refs. iTD, |I2 



3) These RG functions can now be used to calculate the critical exponents describing 
the behaviour of a statistical system near the critical point of the second order phase 
transition ||]. At the critical temperature T = Tq, the asymptotic behaviour of the 
correlation function for |x| — ► 00 has the form 

r(x) ~ . (14) 

|x| 

Close to Tc, the correlation length behaves for t = T — Tq — >• as 

i ~ r"(l + const • f " + ...) . (15) 

The three critical exponents 77, z/ and u defined in this way completely specify the 
critical behaviour of the system. All other exponents can be expressed in terms of 
these [§. 

The three critical exponents can be determined from the RG functions of the 
4 -theory by going to the infrared-stable fixed point 

00 

9 = 9o(e) = Y,9 ik) e k (16) 

k=l 
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which is determined by the condition (/3 e = (3 — |g) 

FM = o, AOto) = m(g)/dg] g=go > . (17) 

The resulting formulas for the critical exponents are: 

?? = 27^0), l/i/ = 2(1 - 7m(pb)) , w = 2pM, (18) 
each emerging as an e-expansion up to order e 5 . From flSp-flTOD we therefore find: 

V® = |w{ 1+ 4i^[-™ 2 + 56 ™ + 272 ] 

- 16(n ^ 8) 4 [5n 4 + 230n 3 - 1124n 2 - 17920n - 46144 

+ C(3)(n + 8) -384(5n + 22)] 
- 64(n £ ^ 8)ti [13n 6 + 946n 5 + 27620n 4 + 121472n 3 

-262528n 2 - 2912768n - 5655552 

- C(3)(n + 8) • 16(n 5 + 10n 4 + 1220n 3 - 1136ra 2 

-68672n - 171264) 
+ C(4)(n + 8) 3 • 1152(5n + 22) 

- C(5)(n + 8) 2 • 5120(2n 2 + 55n + 186)] } , (19) 

l/u(e) = 2 + ^{-1-^(13^ + 44) 
+ 8(n + 8) 4 [3n 3 - 452n 2 - 2672n - 5312 

+ C(3)(n + 8) -96(5n + 22)] 
+ 32 ( n+ 8) e [3™ 5 + 398n 4 - 12900n 3 - 81552n 2 - 219968n - 357120 
+ C(3)(n + 8) • 16(3n 4 - 194n 3 + 148n 2 + 9472n + 19488) 
+ C(4)(n + 8) 3 -288(5n + 22) 

- C(5)(n + 8) 2 • 1280(2n 2 + 55n + 186)] 
+ 128( ^ +8) h [3n 7 - 1198n 6 - 27484n 5 - 1055344n 4 

-5242112n 3 - 5256704n 2 + 6999040n - 626688 

- C(3)(n + 8) • 16(13n 6 - 310n 5 + 19004n 4 + 102400n 3 

-381536n 2 - 2792576n - 4240640) 

- C 2 (3)(n + 8) 2 • 1024(2n 4 + 18n 3 + 981n 2 + 6994n + 11688) 
+ C(4)(n + 8) 3 • 48(3n 4 - 194n 3 + 148n 2 + 9472n + 19488) 

+ C(5)(n + 8) 2 • 256(155n 4 + 3026n 3 + 989n 2 - 66018n - 130608) 
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- C(6)(n + 8) 4 • 6400(2n 2 + 55n + 186) 

+ C(7)(n + 8) 3 • 56448(14n 2 + 189n + 526)] } , 



(20) 



= 6-^[9n + 42] 



+ 4(n + 8) 4 [33n 3 + 538n 2 + 4288n + 9568 

+ C(3)(n + 8) -96(5n + 22)] 
+ 16(n £ + 8)li [5ra 5 - 1488n 4 - 46616n 3 - 419528n 2 - 1750080n - 2599552 

- C(3)(n + 8) • 96(63n 3 + 548n 2 + 1916n + 3872) 
+ C(4)(n + 8) 3 • 288(5n + 22) 

- C(5)(n + 8) 2 • 1920(2n 2 + 55n + 186)] 
+ 64(n ^ g)8 [13n 7 + 7196n 6 + 240328n 5 + 3760776n 4 

+38877056n 3 + 223778048n 2 + 660389888n + 752420864 

- C(3)(n + 8) • 16(9n 6 - 1104n 5 - 11648n 4 - 243864n 3 

-2413248n 2 - 9603328n - 14734080) 

- C 2 (3)(n + 8) 2 • 768(6n 4 + 107n 3 + 1826n 2 + 9008n + 8736) 

- C(4)(n + 8) 3 • 288(63n 3 + 548n 2 + 1916n + 3872) 

+ C(5)(n + 8) 2 • 256(305n 4 + 7386n 3 + 45654n 2 + 143212n + 226992) 

- C(6)(n + 8) 4 • 9600(2n 2 + 55n + 186) 

+ C(7)(n + 8) 3 • 112896(14n 2 + 189n + 526)] . (21) 
For n — 1, these expansions read, numerically, 

r] = 0.01852 e 2 + 0.01869 e 3 - 0.00833 e 4 + 0.02566 e 5 , (22) 

I = 2 - 0.333 e -0.1173 e 2 + 0.1245 e 3 - 0.307 e 4 + 0.951 e 5 , (23) 

to = e - 0.630 e 2 + 1.618 e 3 - 5.24 e 4 + 20.75 e 5 . (24) 

Note that r/^ has decreased by about 30 % in comparison with the (incorrect) results 
of Ref. jnj, has increased by about 10 %, and cu^ increased by about 0.6 % in 
comparison with |T1 . 

It is known that the series of the e-expansion are asymptotic series and special 



resummation techniques |20|, [H[] should be applied to obtain reliable estimates of 
the critical exponents. Although the size of the e 5 terms in the physical dimension 
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(i.e., at e = 1) is very large, their contribution to the exponents in the resummed 
series is very small. This is why even large relative changes of the e 5 coefficients 
turn out not to change the critical exponents P21 within the accuracy of previous 
determinations [ 14 , 23| . 
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